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Abstract. We show that for n > 3 the following equivalence problems are essentially 
the same: the equivalence problem for Lagrangians of order n with one dependent and 
one independent variable considered up to a contact transformation, a multiplication by 
a nonzero constant, and modulo divergence; the equivalence problem for the special class 
of rank 2 distributions associated with underdetermined ODEs z' — f{x, y,y' , . . . , y'"'); 
the equivalence problem for variational ODEs of order 2n. This leads to new results such 
as the fundamental system of invariants for all these problems and the explicit description 
of the maximally symmetric models. The central role in all three equivalence problems is 
played by the geometry of self-dual curves in the projective space of odd dimension up to 
^ , projective transformations via the linearization procedure (along the solutions of ODE 

or abnormal extremals of distributions). More precisely, we show that an object from 
one of the three equivalence problem is maximally symmetric if and only if all curves in 



o 



, projective spaces obtained by the linearization procedure are rational normal curves. 



o 

Q I 1. Introduction: three equivalence problems 

. The main goal of this paper is to estabhsh a tight relationship between the following 

I three local equivalence problems in differential geometry: 

(1) equivalence of variational problems of order > 3 with one dependent and one 
independent variable considered up to a contact transformation, a multiplication 
by a constant, and modulo divergence; 
^> ' (2) equivalence of variational ODEs of order > 6 up to contact transformations. By a 

variational ODE (called also variational with multiplier) we mean an ODE which 
I is contact equivalent to an Euler-Lagrange equations for some Lagrangian. 

(3) equivalence of rank 2 distributions associated with underdetermined ODEs z' = 
/(x,y,y',...,yW), n>3. 

In particular, we shall show that equivalence problems (1), (2), and the equivalence prob- 
lem for the particular class of rank 2 distributions mentioned in item (3) are essentially 
the same. In particular, there is a one-to-one correspondence between equivalence classes 
of objects in all these problems. 

The one-to-one correspondence (up to above equivalence relation) between Lagrangians 



o 



o 
o 



- 

X 



' of order n > 2 and their Euler-Lagrange equations was already established earlier in works 



of M. Fels [16] for n = 2 and M. Juras [20] for n > 3. It is based on the characterization of 
variational ODEs in terms of the variational bicomplex given in [3]. We extend this cor- 
respondence to underdetermined ODEs and the corresponding rank 2 vector distributions 
in the case n > 3. This allows us to apply the results of our previous works [HI US], where 
more general rank 2 distributions are treated, for the description of the unique maximally 
symmetric Lagrangian up to the considered equivalence relation (see also the discussions 
on various equivalence relations for variational problems at the end of subsection II. ip . 
Note that the one-to-one correspondence between the equivalence problems (1) and (3) 
does not hold for n = 2. For example, the Lagrangian {y")^^^ dx is not equivalent to the 
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most symmetric one (y")^ dx, but the corresponding underdetermined ODEs and rank 2 
distributions are equivalent and have 14-dimensional Lie algebra G2 as their symmetry. 

The common feature of all three problems above is that they admit linearization, 
which reduces them in essence to the problem of equivalence of self-dual curves in odd- 
dimensional projective spaces up to projective transformations. The invariants of these 
curves in projective spaces, the Wilczynski invariants [31j . produce the invariants of the 
original problem. The latter are called the generalized Wilczynski invariants. 

In this work we exploit an alternative (a Hamiltonian) point of view on the variational 
problems, which comes from the Pontryagin Maximum Principle in Optimal Control. This 
point of view provides us with the Hamiltonian form of the Euler-Lagrange equation and 
allows to construct a (generalized) Legendre transform that takes extremals of the La- 
grangian (or, in other words, the solutions of the corresponding Euler-Lagrange equation) 
to the abnormal extremals of the corresponding optimal control problem. This immedi- 
ately shows that the solution space of any variational ODE carries a natural symplectic 
structure u. We show that the conformal class of this symplectic structure (i.e., all 2-forms 
fu for non-vanishing functions /) can be recovered only from the self-duality of all lin- 
earizations of the given ODE along its solutions. This, in its turn, can be reformulated by 
vanishing of the generalized Wilczynski invariants (see |31j ) of odd order. It is easy to see 
that for n > 2 there is at most one (up to constant) closed 2-form in any given conformal 
class of non-degenerate 2-forms on a smooth manifold. This gives a 'naive' prove that 
any variational ODE of order > 4 admits at most one Lagrangian up to a constant and 
divergence terms. 

Another question we try to answer in this paper is whether invariants of the above three 
equivalence problems derived from the Wilczynski invariants of self-dual projective curves, 
provide the complete system of fundamental invariants. It has been known from |12[ 
[T3] that the answer to the similar question for arbitrary non-linear ODEs is negative. 
Namely, there exist non-trivial ODEs (i.e., equations, not equivalent to y^"^ = via contact 
transformations) such that all their linearizations are trivial and, thus, all their Wilczynski 
invariants vanish. However, surprisingly, for variational ODEs the answer is positive and 
Wilczynski invariants of even order provide the complete system of fundamental invariants 
for this class of equations. 

We show that for variational ODEs the generalized Wilczynski invariants of even order 
(Wilczynski invariants of odd order automatically vanish due to the self-duality of the lin- 
earization) form a fundamental set of contact invariants in the following sence. Any other 
differential invariant lies in a radical of a differential ideal generated by these invariants. In 
particular, vanishing of this fundamental set of invariants implies that any other differen- 
tial invariant of the variantional equation vanishes identically and the equation is contact 
equivalent to the tivial one. We note that for n = 3 and 4 these invariants do not generate 
the complete differential algebra of invariants by only differentiation and algebraic opera- 
tions. In particular, for n = 3 (or 6-th order variational ODEs y^^^ = F (x , y , y' , . . . ,y^^^)) 
the differential invariant / = F45 satisfies a non-trivial qubic equation, whose coefficients 
are certain derivatives of the generalized Wilczynski invariant of order 4. 

To summarize vanishing of generalized Wilczynski invariants of even order gives an 
explicit characterization of the most symmetric models in all three equivalence problems 
provided that n > 3: 

(1) all variational ODEs with vanishing generalized Wilczynski invariants are contact 
equivalent to the trivial equation y(2n) _ q. 

(2) all Lagrangians with vanishing generalized Wilczynski invariants are equivalent 
to (y("))^(ix modulo constant multiplier, contact transformations and divergence 
terms; 
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(3) all underdetermined ODEs with vanishing generalized Wilczynski invariants equiv- 
alent to z' = (y("))2. 
Let us we briefly outhne each of equivalence problems (l)-(3). 

1.1. Equivalence of variational problems. This paper deals with variational problems 
in one dependent and one independent variable of arbitrary order. A variational problem 
is defined by a Lagrangian L = f(x,y,y',...,y^^^)dx or the corresponding functional 

Let us recall basic definitions from the geometry of variational problems. Let J°° = 
J°°(M, M) be an infinite jets of smooth functions y{x). We shall use the standard coordinate 
system {x,y = yo,yi,y2, ■ ■ ■) on J°°. Denote by 9i = dyi — yi+idx, i > 0, the basis of 
so-called contact forms on J°°. The set of all exterior forms A{J°°) is naturally turned 
into the bi-graded algebra with: 

AP(J°°) = A°'P e A^'P'\ where 

A°'^=(0n A---Aa,^),p>0; 

A^'P-^ = {9,^A... A9,^_^Adx),p>l. 

The exterior derivative d: A^(J°°) — ?> Ap+^(J°°) naturally splits into the sum d = du + dy, 
where for oj G AP(J°°) we have duii^) G A^-^ and dv{oj) E A°'P+^ 

We consider variational problems / f {x , y , y' , . . . , y^""^) dx oi arbitrary order n up to the 
divergence equivalence and constant multiplier. Namely, we say that two Lagrangians 
Li = fi{x, y,y' , . . . , y^"^) dx and L2 = f2{x, y,y' , . . . , y*-"^) dx are equivalent if there exists 
a contact transformation 4>: J^(IR,]R) J^(R,M) with the prolongation <!>: J°°(M,]R) 
J~(M,M), such that 

(1.1) <i>*{L2) = aLi + dnifi) mod (61^ | i > 0) 

for some constant nozero a € M and function on J°°(M,M). We shall always assume that 
all our Lagrangians are non-degenerate, i.e., they are non- linear in the highest derivative. 
It follows from [20] that two Lagrangians are equivalent under the above equivalence 
relation if and only if their Euler-Lagrange equations are contact equivalent. 

The variational equivalence problem was treated in a number of papers using both naive 
approach and Cartan's equivalence method [3 ETJ |22l [231 Ell [13 • See also [27] for the 
symmetry classification of higher order Lagrangians. 

We note that usually slightly different equivalence notion of divergence equivalence is 
considered, where the constant a above is equal to 1 identically. The upper bound for 
the variational symmetry algebra in case of n-order Lagrangian was proved to be equal 
to 2n -|- 3 for n > 2 in the work of Gonzalez-Lopez [IT]. Note that this upper bound is 
achieved in case of a family of non-equivalent Lagrangians, which is different from most of 
the classical local equivalence problems in differential geometry. We show in this paper, 
that adding this constant a in the definition (jl.ip of the divergence equivalence changes 
this patten. In this case we get a slightly higher upper bound equal to 2n + 5 with a 
unique maximally symmetric Lagrangian equivalent to L = {y^"''))'^ dx. 

1.2. Equivalence of rank 2 vector distributions. By a rank 2 vector distribution on a 
smooth manifold M we understand a two-dimensional subbundle D of the tangent vector 
bundle TM. We define its (small) derived flag {D^} as follows: 

= D, 

D'+^ = D' + [D,D% i>l, 

and assume that the distribution D is regular in a sense that all are smooth subbundles 
of the tangent bundle TM. We shall also assume in the sequel that the distribution D is 
completely non-holonomic, i.e. = TM for some sufficiently large n. 
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We say that two such distributions D and D' on manifolds M and M' are (locally) 
equivalent if there exists a (local) diffeomorphism (p: M M', such that (p*{D) = D'. 

Equivalence problem for non-holonomic distributions is an old problem, which goes 
back to the end of 19th century and was studied by various mathematicians including 
Lie, Goursat, Darboux, Engel, Elie Cartan and others. Except for several cases such 
as rank 2 distributions on 3- and 4-dimensional manifolds, generic rank 2 distributions 
have functional, and, thus, non-trivial differential invariants. In his classical paper [B] 
Elie Cartan associates a (2, 5)-distribution to a system of partial differential equations of 
second order and constructs a canonical coframe for non-degenerate distributions of this 
type. This is a first example of an explicit solution for the equivalence problem of vec- 
tor distributions with non-trivial functional invariants. Remarkably, the most symmetric 
(2, 5)-distributions form one equivalence class and have an exceptional Lie algebra G2 as 
their symmetry algebra. 

The obvious (but very rough in the most cases) discrete invariant of a distribution D at 
a point q is the so-called small growth vector (s.g.v.) at q. It is the tuple {dimD^ {q)}j<^fq, 
where . Furthermore, at each point q G M, where dimD^ are locally constant for 
any j, we can consider the graded space rriq = X]j>]^ D^'^^{q) / {q). It can be naturally 
equipped with a structure of a graded nilpotent Lie algebra and it is called a symbol of 
the distribution D at a point q. The notion of symbol is extensively used in works of 
N. Tanaka and his school ([29l [301 [251 [32j) who systematized and generalized the Cartan 
equivalence method. 

However, all constructions of Tanaka theory strongly depend on the algebraic structure 
of the symbol and they were carried out under the very restrictive assumption that symbol 
algebras are isomorphic at different points. An alternative approach for studying rank 
2 distributions was presented by the authors in [lH [15]. It is based on the ideas of 
the geometric control theory and uses a symplectification of the problem by lifting the 
distribution D to the cotangent bundle. This provides an effective way to construct a 
canonical coframe and, thus, solve equivalence problem for rank 2 distributions of so- 
called maximal class (this notion is defined in section [53]) . 

Rank 2 distributions of a special type are naturally associated with Lagrangians. Namely, 
to a variational problem with a Lagrangian f(x, y,y' , . . . , y^"^) dx one can assign the fol- 
lowing (affine) control system: 

x{t) = 1 

^^2^ mit) = yi+i{t), 0<i<n-l 

yn{t) = u{t) 

= f(.x{t),yo{t),yi{t),...,yn{t)), 

on J"(IR, M) xM with coordinates (x, yo, ■ ■ ■ y^ z), where u(-) is a control function belonging 
say to the space L^. To any point q^ G J"'(M,]R) x R and a control function u{-) consider 
the trajectory of the system (jl.2p started at q^. Such trajectory is called an admissible 
trajectory of control system ()1.2p and its velocity at q^ is called an admissible velocity of 
control system ()1.2p at q^. Taking the linear span of all admissible velocities of ()1.2p we 
get the rank 2 distribution on J"(M,M) x M generated by the following two vector fields: 

d d d ^ d 
Xi = — + yi-^ \ hi/n^ 

(1.3) ^ 

"^2 = • 

oyn 

We say that this rank 2 distribution is associated with the Lagrangian f(x, y,y' , . . . , y*-"-*) dx 
or with underdetermined differential equation 

(1.4) z' = f{x,y,y',...,y(^^) 
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on two unknown functions y{x) and z{x). Such underdetermined ODEs and the related 
geometric structures have been extensively studied in fL9[ \T8\ \T0 \ [26 } [3]. 

It is easy to see that if / (i.e. the Lagrangian satisfies the Legendre condition), 

then diuiD^ = 3, dimZ)* = i + 2 for i = 3, . . . ,n + 1. The case, when / is linear with 
respect to Un is special, since in this case we can reduce the corresponding underdetermined 
equation to the equation of the same type, but of lower order. So, in the sequel we shall 
always assume that function / in the right hand of the equation (|1.4p is always non-linear 
with respect to y^"^ . Note that such rank 2 distributions are of maximal class. Finally we 
cite the main result of [m [15] that will be needed in the sequel: 

Theorem 1.1. For any {2, n + 3) -distribution, n > 2, of maximal class there exists a 
canonical frame on a (2n + 5) -dimensional bundle over M . The group of symmetries of 
such distribution is at most (2n + b)- dimensional. Any (2,n + 3) -distribution of maximal 
class with {2n-\- 5) -dimensional group of symmetries is locally equivalent to the distribution, 
associated with the Lagrangian . The algebra of infinitesimal symmetries of 

this distribution is isomorphic to a semidirect sum o/0[(2,R) and (2n + 1)- dimensional 
Heisenberg algebra n2n+i- 

1.3. Equivalence of ordinary differential equations. We shall also consider the equiv- 
alence problem of scalar ordinary differential equations of the form 

(1.5) y(^+i)=F(x,y,y',...,yW). 

Each such equation can be considered as a hypersurface £ in the jet space J^+-'^(M, M). We 
shall always assume that our equations are solved with respect to the highest derivative, 
so that the restriction of the natural projection 7rjv+i,Ar : J^^^(M,]R) — )• J^(R,R) to the 
hypersurface £" is a diffeomorphism. 

Two such equations £ and £' are said to be contact equivalent, if there exists a contact 
transformation J^(IR,M) J^(IR,R) with a prolongation $^+^ to J^+^(R,R) such 
that ^^+\£) = £'. 

The equivalence problem of ordinary differential equations under contact and point 
transformations is yet another classical subject going back to the works of Lie, Tresse, Elie 
Cartan [7], Chern |8j, M. Fels |16) and others. The complete solution for the equivalence 
problem was obtained in |llj based on the Tanaka theory of geometric structures on filtered 
manifolds [23 1321125] • 

Explicit formulas for the basis in the differential algebra of contact invariants of a single 
ODE of arbitrary order were computed by B. Doubrov [12J. In particular, a part of 
these invariants comes from the linearization of the given ODE. In fact, they coincide 
with classical Wilczynski invariants of linear differential equations formally applied to the 
linearization of a non-linear ODE. See |13| for more details. 

In this paper we are mainly interested in a special class of ordinary differential equa- 
tions, consisting of equation which are contact equivalent to Euler-Lagrange equations of 
variational problems. Such equations are usually called variational (with multiplier) and 
have been studied in many papers [20l [T6]. 

From the general result of Anderson and Thompson [5', Theorem 2.6] it is known that 
a scalar ordinary differential equation of order 2n is variational, if and only if there exists 
a closed 2-form: 

n— 1 2n— j— 1 

(1.6) ^ = J2 Yl AA^^j 

i=0 j=i+l 

where An-i^n 7^ and 6i = dyi — yi^idx are contact forms on the jet space J^"(R, R) 
restricted to the equation £ = {y2n = F{x, yo, . . . , y2n-i)}- 
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1.4. Equivalence of curves in projective spaces. Surprisingly, the central role in all 
three equivalence problems is played by the geometry of self-dual curves in the projective 
space of odd dimension up to projective transformations. 

Let 7 C be an arbitrary curve in the projective space. We shall always assume that 
7 is strongly regular, i.e. its flag of osculating spaces does not have any singularities. In 
particular, 7 itself does not lie in any proper linear subspace of the projective space. We 
shall not assume any distinguished parameter on 7, though there is always a distinguished 
family of so-called projective parameters on 7. 

Let t be an arbitrary parameter on 7 and let eo(t) be such curve in that 7(t) = 

Meo(t). Define ej(t) = eQ''(t). Then i-th osculating space "f^^\t) of the curve gamma at a 
point 7(i) is defined as: 

(1.7) j^'Ht) = {eo{t),...,ei{t)), i = 0,...,N. 

It is easy to see that it does not depend on the choice of the parameter t and the curve 
eo(i). The (A'' — l)-st osculating spaces 7^^"^) define the curve in the dual projective space 
P^'*, which is called a dual curve and is denoted by 7*. We shall call a curve 7 self-dual, 
if there exists a projective mapping P^ — )• p^^* that maps 7 to 7* so that any point x G 7 
is mapped to the point in p^'* annihilating the (A^ — l)-st osculating space 7(^~-'^) to 7 
at X . We summarize the properties of self-dual curves in the following proposition. For 
the proofs we refer to the classical book of Wilczynski [31j ) : 

Proposition 1.1. (1) // the mapping P^ — )■ P^>* sending 7 to 7* exists, then it is 
unique, up to a constant nonzero factor. It defines, a unique, up to a constant 
nonzero factor, non- degenerate bilinear form (3 on the vector space MJ^^^. More- 
over, this form is necessarily skew- symmetric if N is odd and symmetric, if N is 
even. 

(2) In case of odd N the curve 7 is self-dual if and only if there exists a non- degenerate 
skew-symmmetric (i.e. symplectic) form on M^+-'^ such that all osculating spaces 
'j(^+^)/'^(^t) are Lagrangian with respect to this form. 

The invariants of projective curves (up to projective transformations) were also de- 
scribed by Wilczynski ^31]. The algebra of all invariants admits a basis of so-called 
fundamental invariants W3, . . . , Wn+i of order 3, . . . , -|- 1 respectively. They can be 
constructed as follows. As above, let eo(t) be a curve in M^+i such that 7 = Reo{t). If 
the curve 7 is strongly linear, then the vectors ej(i) = e^^t), i = 0,... ,N form a so- 
called moving frame along 7. Then the next derivative e'j^{t) = e'^~^^\t) can be uniquely 
expressed as a linear combination of vectors in this frame. In other words, we have a 
well-defined linear homogeneous differential equation on eo{t): 

(1.8) er'"'^=Piv(t)er + ---+Po(t)eo. 

Since eo(t) is defined up to a scale, we can always fix this scaling factor by the condition 
PN^t) = 0. It is easy to see that this defines eo{t) uniquely up to a contact non-zero 
scale. Next, by reparametrizing the curve et(0), i.e., by changing t to t = A(t) we can also 
achieve pwit) = PN-iit) = 0. This fixes a parameter t up to projective reparametrizations 
t = (at + b)/ {ct -\- d). Wilczynski proves that taking linear combinations of the derivatives 
of the remaining coefficients Pi{t), i = 0, . . . ,N — 2we can form (A^ — 2) (relative) invariants 
of the curve 7 under the group of projective transformations: 

Following Wilczynski, we shall say that an invariant Wi has order i, i = 3, . . . , A^ -|- 1. 
Wilczynski proves that any other projective invariant of 7 can be expressed as a function 
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of invariants Wi and their derivatives. He also shows that in case of odd N the curve 7 
is self-dual if and only if all invariants of odd order vanish identically. Note also that all 
Wilczynski invariants of a curve 7 C vanish if and only if 7 is a rational normal curve, 
i.e. it can be represented as t 1-^ [1 : t : . . . : t^] in some homogeneous coordinates of . 

Self-dual curves 7 in odd-dimensional projective spaces appear naturally in the above 
equivalence problems via the linearization procedure (see section H] for more detail). The 
linearization of ODE along a solution assigns a curve in projective space to the solution 
via identification of linear equations with curves in projective space. If the ODE is vari- 
ational, then the corresponding curves in projective space are self-dual. In the case of 
rank 2 distributions it is not immediately clear what is the analog of solutions and what 
is the linearization procedure. This becomes clear if one considers distributions as the 
constraints for a variational problem and use the Pontryagin Maximum Principle: the 
analogs of solutions of ODE's are so-called abnormal extremals of the distribution and 
the linearization of the flow of abnormal extremal leads to the notion of Jacobi curves 
introduced in [35l [HI [15], which essentially are (or generated by) self-dual curves in a 
projective space. In particular, the invariants of these curves define the invariants of the 
original equivalence problems. For example, as shown in [36], the fundamental invariant 
W4 of self-dual curves in RP^ can be identified with the so-called fundamental tensor of 
rank 2 vector distributions in 5 dimensional spaces discovered by E. Cartan [6]. 



2. Variational problems and rank 2 vector distributions 

The aim of this section is to establish the correspondence between variational problems 
of order n > 3 and special rank 2 vector distributions associated with the underdetermined 
ordinary differential equations (|1.4p of order n. 

Lemma 2.1. Let D be the rank 2 distribution associated with to the underdetermined 
z' = f(x, y,y' , . . . , y*-"-*) of order n > 3. Then the space of all infinitesimal symmetries of 
D lying in is one- dimensional (over M) and is generated by the vector field Z = 

Proof. Let us show that sym(L') D is one-dimensional (over M) and is generated by the 
vector field Indeed, it is easy to check that the space is generated by the vector 
fields: 

d d d J. d 

^1 = ^ + yi^ — \ h yn^ h /tt' 

ox dyo oy n-i oz 

X - ^ X - ^ X - ^ X - ^ 
dyn oyn-1 oyn-2 oz 

Let Y = Ylt=i O'i^i be an arbitrary vector field lying in . Then we have 

5 

[Xi,Y] = [Xi,Y,^iXi] modD = 

i=2 

{X,{as) - + {X,{a,) - as)^ - a,^^ + X^-^^ mod D. 

Thus, the condition Y G sym(D) implies that 02 = 03 = 04 = and Xi{a^) = 0. Further, 
we have 

[X2, Y] = X2{a5)j- - ai- mod D. 

dz dyn-2 

Again, the condition Y € sym(D) implies that ai = and ^2(05) = 0. In particular, we 
see that the function 05 is a first integral of the distribution D. But since D is completely 
non-holonomic, 05 should be a constant. This completes the proof of the lemma. □ 
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Corollary 2.1. Let z' = fi{x, y,y', . . . , y^"'^), i = 1,2, be two underdetermined differential 
equations of order n > 3, and let Di, i = 1,2, be the corresponding rank 2 vector distri- 
butions on M"+^. Suppose that distributions Di and D2 are locally equivalent. Then the 
equivalence mapping cj) maps vector field to c-^ for some constant c E M*. 

Let us identify the space M"'^'^ with the direct product of J"(M,M) with the coordinates 
{x,yQ, . . . ,yn) and M with the coordinate z and consider any equivalence mapping as a 
mapping from J"(IR,M) x M to itself. Then Lemma |2 . 1 1 implies that any such mapping (j) 
has the form 

(2.1) (/.: J"(IR,M) X M ^ J"(M,M) X M, (p, z) ^ {i;{p),az + n{p)), 

where ^p: J"(R,M) J"(M,M), a G M* and ^ is a smooth function on J"(M,M). 

Lemma 2.2. Mapping ijj is a contact transformation and the function fi does not depend 
on yn, i.e., it is a pull-back of the function on J"^"^(M,M). 

Proof. Since ^ is a symmetry of both distributions Di and D2, we can consider the 
direct images of these distributions with respect to the natural projection J"(R,M) x M — ?> 
J"(M, M). It is easy to see that in both cases these images coinside with the contact 
distribution on J"'(R,M). This proves that V' is a contact transformation. 

The second statement of the lemma on the function /i follows immediately from the fact 
that both Di and Do contain the vector field 75^. □ 

Theorem 2.1. Suppose that n > 3. Two vector distributions Di and D2 associated with 
Lagrangians Li = fi{x,y,y',...,y^'^^)dx, i = 1,2, for n > 3 are equivalent if and only if 
the Lagrangians Li and L2 are equivalent. 

Proof. It is easy to see that the transformations ()2.ip with mappings if) and /i satisfy- 
ing Lemma 12.21 induce the same equivalence relation on distributions Di and D2 as the 
equivalence relation on Lagrangians Li and L2 given by equation (jl.ip . □ 

As a direct consequence of Theorems 11.11 and 12.11 

Theorem 2.2. The dimension of the group of variational symmetries of Lagrangian of 
order n > 3 does not exceed 2n + 5. The Lagrangian with (2n + 5) -dimensional group of 
variational symmetries is equivalent to the Lagrangian (y*-"^)^. The algebra of infinitesimal 
symmetries of the latter Lagrangian is isomorphic to a semidirect sum o/g[(2,]R) and 
{2n + 1)- dimensional Heisenberg algebra n2n+i- 

From the proof of Lemma 12.21 it follows that rank 2 distribution associated with some 
Lagrangian f{x, y,y', . . . , y*-"-*) dx with ^ can be described in the following coordi- 
nate free way: 

Proposition 2.1. A rank 2 distribution D is associated with a Lagrangian f{x, y, ■ ■ ■ , y*-""-*) dx 
with 0-/0 in a neighborhood of a generic point if and only if 

(1) dimD^ = 5; 

(2) There exists an infinitesimal symmetry X of D lying in such that the factor- 
ization by the foliation of integral curves of X sends D to the Goursat distribution 
on the quotient manifold. 

3. Two POINTS OF VIEW ON EXTREMALS OF VARIATIONAL PROBLEMS 

In this section we introduce abnormal extremals of rank 2 distributions and show 
how the flow of abnormal extremals of a distribution associated with a Lagrangian L = 
f{x, y,y' , . . . , y^"^) dx can be related to the flow of extremals of the corresponding varia- 
tional problems. Speaking informally, this relation is the relation between the Hamiltonian 



EQUIVALENCE OF VARIATIONAL PROBLEMS OF HIGHER ORDER 



9 



and the Lagrangian approach to variational problems and it is given by a kind of Legendre 
transform. The material of this section is pretty standard but, as we shall see in the next 
sections, it is very useful for the equivalence problem for Lagrangians and to our knowledge 
it was never used before in this kind of problems. 

3.1. Hamiltonian form of Euler-Lagrange equation. Recall that extremals of the La- 
grangian L are critical points of the corresponding functional L = J f{x,y,y',...,y^"'^)dx. 
On one hand, they are solutions of the Euler -Lagrange equation 

(3.1) ^^o-^(^^^) + --- + (-l)"^(^^") = 0- 

If one takes a little bit more general point of view (that is standard in the Optimal 
Control Theory), then the extremals can be also described using the notion of the end- 
point mappings associated with the corresponding control system (|1.2p . Fix a point qo G 
J"'(M,M) X M and a time T > 0. The endpoint map J-q^^T is the map from the space 
Loo[0, T] to J"'(M,]R) X M sending a control function u{t) to the point of the corresponding 
admissible trajectory of the system ()1.2p at time T. Then y{t) is an extremal of the 
Lagrangian L if and only if the corresponding control u{t) = y^"'~^^\t) is a critical point 
of the endpoint map J-q^^T for some T > (and therefore for any T > as long the 
corresponding trajectory is defined on [0, T]), where qo = (0, y(0), . . . , y"(0), zq) and zq is 
an arbitrary constant. Take the admissible trajectory q{t) of (|1.'2|) corresponding to the 
control u(t) and starting at qo. Then this trajectory can be lifted to the cotangent bundle 
r*(J"(M,R) X R) by choosing for any t E [0,T] an appropriately normalized covectors 
p{t) G r^*^-|(J"(]R, M) X M) that annihilates the image of the differential dJ^Q^t («(•)) of the 
endpoint map J^qg^t at ^(O- This lifting constitutes one of the main fundamental ideas 
behind the Pontryagin Maximum Principle in Optimal Control ([28j,[lJ). As a matter 
of fact, the curve (p(t),g(t)) is an abnormal extremal of the affine control system ()1.2p 
and also of the distribution associated with the Lagrangian L. This establish in essence 
the relation between extremals of the Lagrangian and the abnormal extremals of the 
corresponding distributions. 

More precisely, the coordinates q = {x,yo, . . . yn, z) in J"(R, R) xR induce the coordinate 
system 

(3.2) {p,q) = (A,^o, • • • ,^n,i^; x,yo,...yn,z) 

in r*(J"(R,R) X R) such that the covector p G r^*(J"(R,R) x R) has the form p = 
Xdx + Y17=o ^ ^ Define the following families of scallar functions (Hamiltonians) 
Hu on r*(J"(R,R) X R): 

n-l 

(3.3) H^{p, q) = \ + Y, + + uf{x, y,y',..., y^")) 

According to the weak form of the Pontryagin Maximum Principle (where the maximality 
condition is replaced by the stationarity condition) on has the following 

Proposition 3.1. A function y{t) is an extremal of the Lagrangian L if and only if for 
the admissible trajectory q{t) of (jl.2p corresponding to the control u{t) = y("+^)(t) and 
starting at the point qo = (0, ?/(0), . . . ,y^"\0),zo), where zq is an arbitrary constant, there 
exists a curve of nonzero covectors p{t) G T*^^^{J"'(K,M.) x R) such that 

(3-4) ^-f^ii(p(i)) 9(i))|M=M(t) = a.e. ^„(t) = (the stationarity condition) 

(3.5) Hm^{p(t),q{t)) = (the transversality condition) 

(3.6) p{t) = —■^H^(i^{p{t),q{t)) (the adjoint equation) 
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Note that another part of the Hamiltonian system 
(3.7) m = ^H^{tMt),q{t)) 

is exactly the system ()1.2p with u{t) = u{t) i.e. it holds automatically. So equations 
()3.4p - ()3.7p can be considered as the Hamiltonian form of the Euler Lagrange equation.. 
The curve {p{t),q{t)) C r*(j"(M,M) X M) satisfying Proposition lO is called an abnor- 
mal extremal of affine control system (jl.2p . The term "abnormal" comes again from the 
Pontryagin Maximum Principle applied to a functional defined on the set of admissible 
trajectories of system (jl.2p : abnormal extremals are exactly the Pontryagin extremals of 
this problem with vanishing Lagrange multiplier near the functional ( |28t [1]). Roughly 
speaking, the extremals of our original variational problem given by the Lagrangian L 
become abnormal extremals of the system ()1.2p . because we include the Lagrangian L 
into this system so that it appears as a part of the constraints. 

Let us analyze the equations ()3.6p in coordinates ()3.2p . First of all, since the Hamiltoni- 
ans ()3.3p do not depend on z we have z> = i.e. v is constant along an abnormal extremal. 
If = 0, then from other equations of (|3.6p and equation ()3.5p it follows that p{t) = but 
p{t) can not vanish by Proposition [3Tj So the case = is impossible. Now assume that 
1/ 7^ 0. From the homogeneity of the equations p.6p with respect to p it follows that it 
is enough to consider the case when u = —1. Then, combining the stationarity condition 
(|3.4p with the equation from (|3.6p regarding ^„ we will get that 



(3-8) C„_l = /y„ 

Writing equations for others from ()3.6p we get 
(3.9) 

ij = fyj - 0-1, J = 1, • • • , - 1; 

Combining ()3.8p and the equation in ()3.9p corresponding j = n — 1 we get S,n-2 = fn-i - 
■^ifn)- Then using the second line of (|3.9p by induction with respect to j in the decreasing 
order, we get 



(3.10) 0-i = E(-l)'"'^(/^'=)' l<i<^-l- 

k=j 

Finally substituting (|3.1Up with j = 1 into the first line of (|3.9p we get the Euler-Lagrange 
equation (|3.ip . as expected. 

If X is a vector field without stationary points or a line disribution, denote by Fol(X) the 
one dimensional foliation of integral curves of X. Consider a codimension 4 submanifold 
n of r*(j"(M,M) X M) given by equations ([231), (1231), and i/ = -1. It is foliated by 

abnormal extremals of system (jl.2p with = — 1. Besides by constructions the group of 
translations along z-axis, z z+c, preserves this foliation. Therefore this foliation induces 
the one-dimensional foliation on the quotient manifold Ti/Fol (^) by the foliation Fol (^) 
or , equivalently, on the manifold of the orbits of the group of these translations. The tuple 
(x, yo, • • • , yn, Co, • • • , Cn-2) constitute a coordinate system on the manifold H/Fol (^) . On 
the other hand, the Euler-Lagrange equation (j3.ip defines a codimension one submanifold 
f (L) of J^"(R, R) foliated by a one- dimensional foliation of prolongations of its solutions 
to J^"'(R, R), and the tuple (x, yo, • • • , 2/2n-i) constitute a coordinate system on £{L). This 
foliation is called the foliation of solutions of the Euler-Lagrange equation. 

By above, the map £: 8{L) ^ n/¥o\ (^), defined by 

(3.11) {x,yQ, . . . ,y2n-i) 1-^ {x,yo, ...,?/„, Co, • • • ,Cn-2), 
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with satisfying ()3.10p , sends the one-dimensional fohation on EL to the one-dimensional 
foliation on "H/Fol (^)- In other words, this map transforms the extremals of our vari- 
ational problem obtained in the Lagrangian form to the extremal obtained in the Hamil- 
tonian form. Therefore we call it the (generalized) Legendre transform. Note that the 
Legendre transform depend on the choice of coordinates on J'^(M,M) = M^, which induces 
the coordinates on J^"(M, M)). Once we use the Legendre transform in the sequel it will 
mean that such choice is already done. 

3.2. Abnormal extremals of rank 2 distributions. Now we are going to describe 
abnormal extremals for a distribution D on a manifold M. We shall use more geometric 
language. Let vr : T*M i— )• M be the canonical projection. For any A S T*M, X = {p,q), 
q G M, p S T*M, let s(A)(-) = p{Tr^-) be the tautological Liouville 1-form and a = ds he 
the standard symplectic structure on T*M. Denote by ^ T*M the annihilator of 

the Zth power D\ namely 

(3.12) {D^)^ = {{q,p) G T*M : p • u = V?; G D\q)}. 

Finally let Sq be the zero section of T*M. With this notation the Pontryagin Maximum 
Principle in the coordinate- free form ([1]) implies immediately the following description of 
abnormal extremals of the distribution D: 

Definition 3.1. An absolutely continuous curve 7 C T*M is an abnormal extremal of a 
distribution D if the following two conditions holds: 

L 7 c D^\So, 

2. 7(f) belongs to Ker(cj|^x); ^-c, to the kernel of the restriction of the canonical 
symplectic form a to the annihilator of D. 

From now on we will consider only rank 2-distributions. From direct computations |33| 
Proposition 2.2] it follows that Ker(c7(A)|^x) 7^ if and only if A G (-D^)-*-. This implies 
the following characterization of abnormal extremals of rank 2 distribution. 

Proposition 3.2. An absolutely continuous curve 7 C T*M is abnormal extremal of a 
rank 2 distribution D with dimD^ = 'i if o-nd only if the following two conditions holds 

L r C {D^)^\So, 

2. T{t) belongs to Ker((j|^^2)±); ^-e., to the kernel of the restriction of the canonical 
symplectic form a to the annihilator (-D^)"*" of . 
Further, if A G {D'^)-^\{D'^)-^ then Ker(o"|^j-,2')±) is one-dimensional. These kernels form 

a special line distribution on A G {D'^)^\{D'^)^ , which will be denoted by C, called the 
characteristic distribution. The abnormal extremals of D, lying in A G {D'^)-^\{D'^)^ , are 
exactly the integral curves of the line distribution C (in some literature these abnormal 
extremals are called regular). 

Remark 3.1. For any A G {D'^)^\{D^)^ let 

(3.13) J{\) = {ve Tx{D^)^ : TT.ve D{Tr{X)) }. 

A simple count shows that dim J7(A) = n + 2. Then from constructions it follows imme- 
diately that the restriction of the form cr(A) to ^^(A) is identically equal to zero. 

Now consider the distribution D on M = J"(M,M) x M associated with the Lagrangian 
L = f{x,y,y', . . . with 7^ 0. We would like to rewrite the constructions of 

the end of the previous subsection in more geometric form. First of all in the considered 
case Z)-*" is a corank 2 submanifold of T*M given by equations (|3.4p and (|3.5p . (-D^)-*- 
is a corank 1 submanifold of satisfying in additional equation (j3.8p . and (-D^)"*" is a 
corank 2 subdistribution of {D"^)-^ satisfying to additional equations u = and Cn-2 = 0. 
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The submanifold T-L introduced in the previous subsection is equal to {Hz = —1} H (-D^)"*" 
and the abnormal extremals of the distribution D lying on 7i coincide (as unparametrized 
curves) with the abnormal extremals of system ()1.2p having = — 1. 

Further, given a vector field X on M denote by Hx : T*M ^ M. the corresponding 
quasi-impulse 

Hxip, q) = p{X{q)) , g e M,p G T^M, 

and by Hx the corresponding Hamiltonian vector field on T*M, i.e. the vector field satis- 
fying ifj^cr = —dHx- It is clear that if X is an infinitesimal symmetry of the distribution 

D, then the flow e*^^ , generated by Hx^ sends an abnormal extremal of D to an abnormal 
extremal of D. Moreover, any abnormal extremal lies on a level set of the function Hx- 
In particular, let as in Lemma [2. II Z be the inflnitesimal symmetry of D lying in D^. The 
distribution C induces a rank 1 distribution C on the quotient manifold T-L/Fol{Hz), where 
as before Fol{Hz) is the foliation of integral curves of the field Hz- A Legendre transform 
£: £{L) H/ Fol{Hz), defined in the previous subsection, sends the one-dimensional 
foliation of solutions of Euler-Lagrange equations to the one-dimensional foliation of the 
integral curves of the distribution C- 

Remark 3.2. The Legendre transform £ satisfies another important property. To describe 
it in geometric terms let iTij : J*(R,M) — > J-'(IR, M) , where i > j, and 

7f : r*(J"(M,M) X M) ^ J"(M,M) x M 

denote the canonical projections. The mapping vf induces the mapping 

1TZ-- T*(J"(M,M) X R)/Fol{Hz) ^ (J"(M,M) x M)/Fol(Z) ~ J"(IR,M) 

in the obvious way. If the submanifolds £{L) and 7i/Fol{Hz) are considered as fiber 
bundles over J"(]R, M) with the projections iT2n,n\i:i^L) ^^l-H/e*^z' ^'^spectively, then 
from (|3.1ip it follows immediately that the Legendre transform £ is fiberwise mapping 
over the identity on the base manifold J"(R,M). 

Remark 3.3. The tautological Liouville 1-form s and the standard symplectic structure a 
on T*(J"(M,M) X M) induce the 1-form s and the closed 2-form a = ds, respectively, on 
Ti/ Fol{Hz) ■ By constructions, the rank 1 distribution C satisfies C = Keia. Besides, 
using condition ()3.5p it is easy to show that 

n-2 

(3.14) s = L + Y,^^Si + fyJn-l 

i=l 

in the coordinates (x, yo, . . . , yn,Co, ■ ■ ■ , Cn-2) on H/ Fol{Hz), where, as before, 6i = dyi — 
yi+i dx. Finally, for any A € 'H/Fo\{Hz) we denote 

(3.15) J{\) = {ve Txn/Fol{Hz) : {7Tz)*v (£ D{tt{X))}. 

Then from the last sentence of Remark l3. li the restriction of the form (t{\) to the subspace 
j7(A) is identically equal to zero. 

4. Linearization of variational ODEs and Jacobi curves of rank 2 

distribution 

Let us outline the content of this section. As was mentioned before both in the equiva- 
lence problem for variational ODEs and in the equivalence problem for rank 2 distributions 
self-dual curves in a projective space play a crucial role. They appear via the linearization 
along the "flow "of solutions in the flrst case and along the "flow" of abnormal extremals in 
the second case. Using the Legendre transform introduced in subsection 13 . 1 1 we show that 
the self-dual curves in a projective space obtained by the linearization along a solution of 
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the Euler-Lagrange equations of a Lagrangian and by the hnearization along the corre- 
sponding abnormal extremal of the associated rank 2 distribution are actually isomorphic. 
This observation leads to the description of the fundamental system of invariants for rank 
2 distributions associated with Lagrangians given in the next section. 

4.1. General linearization procedure. Let us first clarify what do we mean by the 
linearization procedure in a general geometric setting. Let M be an arbitrary smooth 
manifold, let G and V be a pair of vector distributions on Ai of rank / and k, respectively, 
where one of them, say G, is integrable and V CiG is a distribution of rank r. 

Similarly to above, let Fol(^) be a foliation of A4 by maximal integral submanifolds of 
G- Then we can define the linearization of the distribution V along the foliation Fol{G) 
in the following way. Let, as above, Fol(^) be a foliation of A4 by maximal integral 
submanifolds of G- Locally we can assume that there exists a quotient manifold A4 / Fol(^), 
whose points are leaves of J~{G)- Let F by any such leaf. Then we define the map (f> 
of F into the Grassmannian Grk-riTr{-M./ Fol{G)) of (A; — r)-dimensional subspaces of 
(A1/Fo1(^)) or, under additional regularity assumptions, an /-dimensional submanifold 
of Grfc_^(rr(M/Fol(^))) as follows: (j){x) = pr,(V^), x G F, where pr: M M/Fol{G) 
is a natural projection. The map or its image in Grk-r(Tr(A^/Fol(t/))) is called the 
linearization of the distribution V along the foliation Fol(^) at the leafT or the linearization 
of the distribution V along the leafT (of Fol{G))- In the cases under consideration / = 1 
so that the linearizations are curves in projective spaces. 

The main idea of using the linearization procedure in the equivalence problem for the 
structures given by the pair of distribution (V, G) on A4 with respect to the action of group 
of diffeomorphisms of M is that it allows to construct the invariants of such structures 
from invariants of submanifold in an appropriate Grassmannian with respect to the natural 
action of the General Linear Group on this Grassmannian. 

4.2. Linearization procedure for ODEs. As in subsection II. 31 an ODE of order + 1, 
resolved with respect to the highest derivative, is given by a hypersurface £ in the jet 
space J^+-^(R,M) so that the restriction of the natural projection ttn+i,n- 

J^(M,M) to the hypersurface £" is a diffeomorphism. Further, the Cartan distribution of 
/^"•"^(MjR) (i.e. the rank 2 distribution defined by contact forms 9i = dyi — yiJ^idx, i = 
0,...,N, in the standard coordinates (x, yo; • • • ; ^Af+i) in J^'''^(M,R)) defines the line 
distribution iS on This distribution is obtained by the intersection of the Cartan dis- 
tribution with the tangent space to £ at every point of £. Note that the corresponding 
foliation Fol(S') is the foliation of solutions of our ODE (more precisely, the foliation of the 
prolongations of the solutions to J^(IR,M)). If the hypersurface £ has the form ()1.5p in 
coordinates {x,yo, . . . ,yAr+i) on J^^"-'^(R, M), then in the coordinates {x,yo, . . . ,yN) on £: 

(4.1) S=l^- + Y. y^+i 9^ + ^(^' yo, 2/1, . . . , 2/^)^^ . 

The distribution S will play the role of the distribution G from the previous subsection. 

Further let, as before, 7rAr+i_Ar_j : J^^^ (M, M) i— )• J^~*(]R, M) be the canonical projection. 
For any e € £" we can define the filtration {V^*}^q of T^£ as follows: 

(4.2) = ker deTTN+i,N-i n Te£. 

Then is a rank i distribution on In the coordinates {x,yo, . . . ,yj\f) on £ we have 

(4.3) = / 

\ dyN-i+i oyN , 

Let Sol denote the quotient manifold £/Fol{S), i.e. the manifold of solutions of the 
equation £. Fix a point F € Sol. In other words, F is a leaf of Fol(5) or a solution of the 
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equation £. Consider the Unearization Liup of the distribution along F. It is a curve 
in Grj(rrSol). In particular, Liup is a curve in the projective space P(rrSol). Moreover, 
if r is considered as the leaf of Fol(5), then from (|4.ip and (j4.2p one gets immediately the 
following 

Lemma 4.1. For any e G F the i- dimensional subspace Linp(e) 0/ TpSol is exactly the 

1- th osculating space of the curve Lin^ at e (as defined in (jl.7p ). 

The Wilczynski invariants of the linearizations Liup taken for every solution F define 
the invariants of the original ODE £ under the group of contact transformations (see [12j). 
We call these invariants the generalized Wilczinski invariants of £ and denote them also 
hy Wi,i = 3,..., N + 1. 

4.3. The case of variational ODEs. Now assume that £ is a variational ODE, £ = £{L) 
for some Lagrangian L. Let a be the closed 2-form on H/ Yo\{Hz) introduced in Remark 
13.31 and let £: £{L) — > 7i/Fol{Hz) be a Legendre transform. Then we can define the 
closed 2-form lo on £{L) as follows: 

(4.4) u = 2.*a 

Note that from relation (|4.3p and Remark l3.2l it follows that the distribution defined 
by (|3.15p satisfies 

(4.5) J = £*(K)ec. 

From this and Remark 13.31 it follows that the form oj satisfies the following two properties: 

(1) Keioj = S, where, as before, S is the rank 1 distribution generating the foliation 
of solutions of £{L); 

(2) The restriction of uj to the distribution vanishes. 

From property (1) it follows that uj induces the symplectic form a) on the manifolds Sol 
of solutions of £{L). Moreover from the property (2) it follows that 

(2') For any F G £{L) the linearization Linp of the distribution along G is the curve 
of Lagrangian subspaces with respect to the symplectic form w(F). 
From item (2) of Proposition 11.11 we get immediately the following 

Corollary 4.1. For a variational ODE the linearizations Liup along any solution V is a 
self-dual curve in the corresponding projective space. 

Further, as an immediate consequence of item (1) of Proposition 11.11 and the fact that 
UJ is closed, we get 

Proposition 4.1. For a variational ODE £{L) there exists a unique, up to a constant 
nonzero factor, closed 2-form u on £{L) satisfying conditions (1) and (2) above or, equiv- 
alently, a unique, up to a constant nonzero factor, symplectic structure a) on the manifold 
of solutions Sol, satisfying condition (2') above . 

Remark 4.1. It is easy to see that this symplectic strucutre u) from Proposition 14.11 is 
exactly the 2-form oj given by ()1.6p and prescribed by Theorem 2.6 of Indeed, the 
condition duj = Q implies duu) = 0. In particular, this means that oj projects to the 
solution space of the equation £{L). Proposition gives an alternative construction for this 

2- form based on the generalized Legendge transform. 

Now we are ready to prove the following Theorem (which is also proved in \16\ Theorem 
5.1] in the case n = 2 and in [201 Corollary 2.6] for n > 3): 

Theorem 4.1. Lagrangians are equivalent if and only if their Euler-Lagrange equations 
are contact equivalent. 
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Proof. If Lagrangians Li and L2 are equivalent then directly from (jl.ip it follows that 
their Euler-Larange equations are contact equivalent. 

In the other direction let ip be the mapping establishing the equivalence of two Euler- 
Lagrange equations £1 and £2 ■ Let wi and UJ2 be the closed 2-forms from Proposition I4.H 
corresponding to £1 and £2, respectively. Then by this Proposition there exists a constant 
a / such that 

(4.6) ip*uj2 = au)i. 

Now assume that Si and S2 are 1-forms from Remark 13.31 corresponding to Lagrangians 
Li and L2, respectively. Assume that £i are the corresponding Legendre transforms and 
let Pi = (£j)*Sj, i = 1,2. Then by construction oji = dpi. Hence from ()4.6p there exists a 
function p on £1 such that 

■0*P2 = api + dp. 

Taking into account the coordinate expressions for the forms Sj given by (j3.14p . we get 
immediately that the last relation is equivalent to (jl.ip . i.e. the Lagrangians Li and L2 
are equivalent. □ 

4.4. Linearization procedure for rank 2 distributions: Jacobi curves. The presen- 
tation of this subsection is rather closed to our previous works jl4[ 115] but it is considered 
here in relation with the linearization procedure of ODE's from subsection ()4.2p . Let D 
be an arbitrary rank 2 vector distribution on an [n + 3)-dimensional manifold M. We 
shall assume that D is completely non-holonomic and that dimL*^ = 5. Here it is more 
convenient to work with the projectivization of Pr*M rather than with T*M. Here ¥T*M 
is the fiber bundle over M with the fibers that are the projectivizations of the fibers of 
T*M. The canonical projection H: T*M ^ ¥T*M sends the characteristic distribution C 
of {D'^)^\{D^)^ to the line distribution C on F{D^)^\F{D^)^ , which will be also called 
the characteristic distribution of the latter manifold. The manifold P(L'^)^\P(D^)^ and 
the distribution C play the role of A4 and G, respectively, from the general linearization 
procedure of subsection 14.11 

Further note that the corank 1 distribution on T*M\Sq annihilating the tautological 
Liouville form s on T*M induces a contact distribution on PT*M, which in turns induces 
the even-contact (quasi-contact) distribution A on ¥{D'^)^\¥{D^)^ . The characteristic 
line distribution C is exactly the Cauchy characteristic distribution of A, i.e. it is the 
maximal subdistribution of A such that 

(4.7) [C, A] C A. 

Now let J be as in ()3.13p . Let J{X) = Ti^fJ and define a sequence of subspaces J^''-\X), 
A G P(L'^)^\P(D^)^, by the following recursive formulas: 

jW(A) = j{^-i)(A) + [C,^^-i](A). 

By |35^ Proposition 3.1], we have dim^(^^(A) — dim J'(A) = 1, which implies easily that 

dim J^W(A) - dim j(*-^)(A) < 1, Vi G N. 

Besides, from ()4.7p it follows that J^"^^ C A for all natural i. Simple counting of dimensions 
implies that rank A = 2n + 1 so that dim J'^*^(A) < 2n + 1. 

Note that for any A € F{D'^)^\F{D^)^ the subspace A(A) is equipped canonically, up 
to a constant nonzero factor, with a skew-symmetric form with the kernel equal to C(A): 
for this take the restriction to A(A) of the differential of any 1-form annihilating A. Given 
a subspace W of A (A) denote by the skew-symmetric complement of W with respect 
to this form. Note that by RemarkO (^^°^)^ = Then set 

(4.8) J(-*)(A) = (^(-'^(A))^. 
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The sequence of subspaces {J'^''\X)}i^z defines the filtration of A(A). 
Further, define the following two integer-valued functions: 

u{\) = mm{i G N : j'^'+^^X) = j'''\X)}, 

m{q) = max{u{X) : X G {D^)^{q)\{D^)^{q)}, q G M. 

The number m{q) is called the class of distribution D at the point q. By above, 1 < 
m{q) < n. It is easy to show that germs of (2, n + 2)) -distributions of the maximal class n 
are generic ( see [351 Proposition 3.4]). 

From now on we assume that D is a (2, n + 3)-distribution of maximal constant class 
m = n. Let 7^ = {A G F{D'^)-^\¥{D^)-^u{X) = n}. Then on A G 71 the subspaces J^'^ 
form a distribution of rank (i + n + 1) for all integer i between — n and n (in particular 
ji^) = A). If we denote by Abn the quotient of manifold ¥{D^)^\F{D^)^ by Fol(C) 
then the distribution A induces the distribution A on Abn equipped with the canonical, 
up to a constant nonzero factor, symplectic form. Given any segment T of abnormal 
extremal ( a leaf of Fol(C)) consider the linearization J^*^ of the distribution J'^^^ along 
T. It is a curve in Grn+i(A(T)). In particular, is a curve in the projective space 

P(A(r)). The curve Jp is called the Jacobi curve along the abnormal extremal F. By 
Remark 13.11 it is the curve of Lagrangian subspaces of A. Further, it is not hard to see 
([35j) that for any A G F the (n + i)-dimensional subspace J^\X) of TxAbn is exactly the 
i + n — 1-st osculating space of the curve at A (as defined in (jl.7p ). So, by item (2) 

of Proposition 11.11 the curve is self-dual. 

Remark 4.2. It can be shown [34j that is the only curve in P(A(T)) such that the 

Jacobi curve at A G T is its (n — l)-st osculating space at A. □ 

The Wilczynski invariants of the linearizations taken for every abnormal T de- 

fine the invariants of the distribution D. The latter invariants are called the generalized 
Wilczinski invariants of 8. The Jacobi curve along the abnormal extremal is called flat 
if the corresponding curve is a rational normal curve in P(A(T)) or, equvalently, 

all Wilczinsky invariants of vanishes identically. For the maximally symmetric 

(2, n -|- 3)-distribution of maximal class which is locally equivalent to the distribution 
associated with the Lagrangian {y^^^)'^ dx (Theorem ll.ip all Jacobi curves are flat or, 
equivalently, all generalized Wilczynski invariants vanish. The general question is 

Question Is it true that if all generalized Wilczynski invariants o/(2,n + 3) distribution 
of maximal class vanish or , equivalently, all its Jacobi curves are flat, then the distribution 
is locally equivalent to the distribution associated with the Lagrangian {y^"''^)'^ dx? 

Since, as was shown in [36j, in the case n = 2 the generalized Wilczynski invariant 
coincides with Cartan's covariant binary biquadratic form introduced in [6], then the fact 
that this form is the fundamental invariant of a (2, 5)-distribution (proved in [6] as well) 
gives the positive answer to our question in this case. We show in the next section that 
for n > 3 the answer to this question is positive if we restrict ourselves to distributions 
associated with Lagrangians (see Theorem 15.1 1 below) . Note that this class of distributions 
is rather restrictive and the general question for n > 3 remains open. 

Assume that the distribution D is associated with some Lagrangian L with ^ / 0. 
Then in the notations of section the manifold % can be identified with an open subset of TZ. 
Let (j)z ■ T-Lz T~Lz / ^o[{Hz) be the canonical projection. Then, as was already mentioned 
in subsection 13.21 the characteristic distribution C is reduced to the line distribution C = 
{4'z)*C on H/Fol{Hz)- Moreover, the filtration {J^^^^}i^z on 7i induces the filtration 
on H/ Fol{Hz)- Note that in this notation the distribution J'^'^^ coincides with 
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J defined by ()3.15p . If Ti and T2 are two abnormal extremals on H. such that ^z^\) = 
(pz{'^2)- then the curves J^^ and J^^ coincide up to a projective transformation. As 
a matter of fact they coincide, up to a projective transformation, with the linearization 
of the distribution ^(^~"-) along T = {<j)z)*'^i on %/ ¥o\{Hz). From this. Lemma [4.11 
Remark 14.21 ci^id ()4.5p we have the following 

Proposition 4.2. Given a solution T of the Euler- Lagrange equation £{L) and the ab- 
normal extremal T such that (pzC^) = ^i^), where £, is a Legendre transform, the curves 
Liup and coincide up to a projective transformation. 



5. Fundamental invariants for our equivalence problems 
Now we are ready to prove the following 

Theorem 5.1. For any n > 3 a {2,n + 3) -distribution D associated with a Lagrangian 
f{x, y,y', . . . , y^"^) dx with |^ ^ (or the Lagrangian L itself) is locally equivalent to the 

distribution associated with the Lagrangian (y^^^)"^ dx if and only if one of the following 
two equivalent conditions is satisfied: 

(1) all Jacobi curves of the distribution D are flat; 

(2) generalized Wilczynski invariant of D vanish identically. 

The proof of the theorem immediately follows from Theorem 14.11 and the following 

Theorem 5.2. Let y^^"-) = F{x,y,y' , . . . ,y^'^^~^^) be the Euler- Lagrange equation of the 
Lagrangian L = f{x,y,y',...,y^"^)dx. This equation is contact equivalent to the trivial 
equation y^^") = if and only if all its generalized Wilczynski invariants W4, Wq, . . . , W2n 
vanish identically. 

Proof. In the sequel we will denote by Fi the partial derivative Fy.. The higher order 
derivatives of F will be denoted in a similar way. According to [12] , any ordinary differen- 
tial equation of order 2n, n > 3, is trivializable if and only if all its generalized Wilczynski 
invariants vanish identically and, in addition, the following conditions hold: 

(5.1) F55 = F45 = for n = 3; 

(5.2) F2n^l,2n-l = -^2n-l,2n-2 = -^2n-2,2n-2 = forn > 4. 

Since our equation is variational, its generalized Wilczynski invariants of odd degree vanish 
automatically. The generalized Wilczynski invariants of even degree vanish by assumption 
of the lemma. So, we need to prove that the above conditions also hold. 

Since the equation = F{x, y,y', . . . , y(2n-i)~) variational, then according to [3] F 
is a polynomial in y("+^), . . . , of weighted degree < n, where these derivatives have 

weights 1, . . . , n — 1 respectively. In particular, we see that the polynomials (y^^"~^^)^, 
y{2n-i) y{2n~2) ^ (y(2n-2)^2 -j^^^^ weighted degree 2n — 2, 2n — 3 and 2n — 4 respectively. 
Assume that n > 5. Then 2n — 4 > n, and these terms can not appear in F. Thus, the 
condition (j5.2p holds automatically for n > 5. So, it remains to consider only the cases 
n = 3,4. 

Let n = 3. Then the term y^^^ has weighted degree 2, while the function / is of weighted 
degree < 3. So, we see that F55 = 0. Let us prove that the condition VF4 = implies also 
that F45 = 0. The direct computation shows that / = F45 = —3/3337/33. Let us denote 
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W4 simply by W. From |12) we have: 

5 2 5 1 

6b 21 12 6 

+ TZ^^F^xx + i^^F^F^ — TT^F^F^x — -T^T^FiF^x 
18 6b 21 12b 

5 37 25 5 5 

+ 25^^' + 630^' + 84^'- + 18^^ " Ts^^^^^' 

where Fi denotes the partial derivative by y^*-* and Fx denotes the total derivative. Then 
the direct computation shows that: 



W445 



_ 1 57/I33 - 35/33/3333 

35 /I3 



1 35/33/33333 — 149/33/333/3333 + 144/^ 



333 . 



35 /3 



33 

2 



2 35/33/33333 - I62/33/333/3333 + 135/333 
'35 /I 



'33 

and we have the following syzygy: 

2101^355 - W5Wm5 + 26IW55 - T^^^ = 0. 

105 

In particular, if W = we get 1^ = and hence / = 0. In particular, the function / is 
actually quadratic in the highest derivative. This proves the case n = 3. 

For n = 4, we see that F77 = Fjq = due to the weighted degree argument. So, 
it remains to prove that / = Fqq = — 6/4447/44 vanishes if all generalized Wilczynski 
invariants vanish identically. Again, denote hy W = W4, the first non-trivial Wilczynski 
invariant for Euler-Lagrange equation. Then we have: 

1127 2z:, 161 „ ^2 931 2 7 

— Fn -Tyr H Fr F'j H Fyrr- H Fa 

6336 ^ 3168 ^ 3168 66 

47 2 1127 -^4 329 49 

H Fc H F'J — FaFjx — i*7 7.9-7.. 

1584 ^ 101376 ^ 1584 ^ ' 264 ' 

Again, direct computation proves shows that: 

TJ7 7 ?>f\if\m — '^'ifliA 

66 /44 



2 ' 

2 



pr^ 1 252/44/4444 — 437/444 

198 /I4 

and we have the following syzygy: 

3W^75 - 21^66 + tJ^/' = 0. 

648 

Hence, the equality W = Q implies also that / = 0. Again, we see that vanishing of the 
first non-trivial generalized Wilczynski invariant implies that the function / is actually 
quadratic in the highest derivative. This completes the case n = 4 and the proof of the 
Lemma. □ 

Remark 5.1. Note that for n = 2 Theorem 15 . 21 does not hold. In particular, the Lagrangian 
L = (y")^^^ dx has trivial Wilczynski invariants, but the associated Euler-Lagrange equa- 
tion 3y"y^^^ — 5{y"')^ = is not trivializable and has only 6-dimensional symmetry algebra. 
On the other hand. Theorem 15. II is valid for n = 2 and not only for distributions associated 
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with second order Lagrangians but for any rank 2 distribution in with the small growth 
vector (2,3,5). The rank 2 distributions D corresponding to the equations z' = {y"Y^'^ 
and z' = (y")^ are equivalent and have the 14-dimensional symmetry algebra, which is the 
maximal possible algebra for distributions under the consideration. 

Remark 5.2. Direct analysis of the symmetry classification of all ordinary differential 
equations (see [27J) shows that all Euler-Lagrange equations of order 2n with the symmetry 
algebra of dimension at least 2n + 1 are exhausted (modulo contact transformations) by 
the following ones: 



Equation 8 


Lagrangian L 


dim sym(£^) 


y2n = 


y'ridx 


2n + 4 


yl + T,7=o (^iv^i = 


{yl + J:t=oc^y^) dx 


2n + 2 


Syfye - ^^ywm + 40y| = 


^3^^ dx 


7 



Note also that we always have dimsym(L) = dimsym(£^) + 1. 
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